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Abstract 



There is a conformal equivalence between power law f{R) theories and scalar field 
theories in which the scalar degree of freedom evolves under the action of an exponential 
potential function. In the scalar field representation there is a strong coupling of the scalar 
field with the matter sector due to the conformal transformation. We use chameleon 
mechanism to implement constraints on the potential function of the scalar field in order 
bi)! that the resulting model be in accord with Solar System experiments. Investigation of 

these constraints reveals that there may be no possibility to distinguish between a power 

law f{R) function and the usual Einstein-Hilbert Lagrangian density. 
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There are strong observational evidences that the expansion of the universe is accelerating. 
K^ ! These observations are based on type la supernova [1], cosmic microwave background radiation 
H I [2], large scale structure formation [3], weak lensing [4], etc. The standard explanation invokes 
an unknown component, usually referred to as dark energy. It contributes to energy density 
of the universe with Qd = 0.7 where Qd is the corresponding density parameter, see e.g., [5] 
and references therein. The simplest dark energy scenario which seems to be both natural 
and consistent with observations is the ACDM model in which dark energy is identified as a 
cosmological constant [5] [6] [7]. However, in order to avoid theoretical problems [6], other 
scenarios have been investigated. Among these scenarios, there are quintessence [8], tachyons 
[9], phantom [10], quintom [11] and modified gravity models [12]. In the latter, one modifies 
the laws of gravity whereby a late time acceleration is produced without recourse to a dark 
energy component. One family of these modified gravity models is obtained by replacing the 
Ricci scalar R in the usual Einstein-Hilbert Lagrangian density for some function f{R). How- 
ever, changing gravity Lagrangian have consequences not only in cosmological scales but also 
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in galactic ones so that it seems to be necessary to investigate the low energy limit of such 
f{R) theories. 

Early works on the weak field limit of f{R) theories led to negative results. In fact, using 
the equivalence of f{R) and scalar-tensor theories [13] [14] [15], it is originally suggested that 
all f{R) theories should be ruled out [17] since they violate the weak field constraints coming 
from Solar System experiments. This claim was based on the fact that f{R) theories (in the 
metric formalism) are equivalent to Brans-Dicke theory with u = while observations set the 
constraint u > 40000 [18]. In this case the post-Newtonian parameter satisfies 7 = ^ instead 
of being equal to unity as required by observations. Later, it was noted by many authors that 
for scalar fields with sufficiently large mass it is possible to drive 7 close to unity even for null 
Brans-Dicke parameter. In this case the scalar field becomes short-ranged and has no effect 
at Solar System scales. Recently, it is shown that there exists an important possibility that 
the effective mass of the scalar field be scale dependent [19]. In this chameleon mechanism, 
the scalar field may acquire a large effective mass in Solar System scale so that it hides local 
experiments while at cosmological scales it is effectively light and may provide an appropriate 
cosmological behavior. 

In the present work we intend to use this criterion to set local gravity constraints on f{R) 
theories. There are a number of works concerning these constraints on f{R) theories [20] [21]. 
We will focus on power law f{R) theories and show that the constraints on the parameters 
space suggest that they are hardly distinguishable from ACDM scenario. 

Let us begin with the following action''" 

1 r -4. 



S=-J d'^xy^ f\R) + Srn{g^.u. ^) (1) 

where g is the determinant of g^v, f{R) is an unknown function of the scalar curvature R and 
Sm is the matter action depending on the metric g^u and some matter field ip. We may use a 
new set of variables 

gtiu = P g^v (2) 

</'=^lnp (3) 

where p = -^ = f (R) and (3 = w|. This is indeed a conformal transformation which 
transforms the above action in the Jordan frame to the Einstein frame [13] [14] [15] 

S=^J d^x./^ {R - r^^d^cj) 9,0 - 2\/(0)} + S^isi^^e^^'^, ^j) (4) 

In the Einstein frame, is a scalar field with a self-interacting potential which is given by 

1^(0) = \e-^^'^{T[pm - e-^^^/(r[p(0)])} (5) 



■fWe use the unit (SttG)-! = 1. 



where r(p) is a solution of the equation / [r(p)] — p = [13]. One usually states that any f{R) 
gravity model is mathematically equivalent with a minimally coupled scalar field theory with an 
appropriate potential function. In general, this does not imply physical equivalence of the two 
conformal frames. In fact it is shown that some physical systems can be differently interpreted 
in different conformal frames [13] [16]. The physical status of the two conformal frames is an 
open question which we are not going to address here. However we assume that the scalar 
degree of freedom in the Einstein frame should satisfy stringent constraints from solar system 
experiments. It is important to note that conformal transformation induces the coupling of 
the scalar field with the matter sector. The strength of this coupling (3, is fixed to be w | and 
is the same for all types of matter fields. In the case of such a strong matter coupling, the role 
of the potential of the scalar field is important for consistency with local gravity experiments. 
When the potential satisfies certain conditions it is possible to attribute an effective mass to 
the scalar field which has a strong dependence on ambient density of matter. A theory in 
which such a dependence is realized is said to be a chameleon theory [19]. In such a theory 
the scalar field can be heavy enough in the environment of the laboratory tests so that the 
local gravity constraints suppressed even if /3 is of the order of unity. Meanwhile, it can be 
light enough in the low-density cosmological environment to be considered as a candidate for 
dark energy. 
Variation of the action (1) with respect to g^j^iy and 0, gives the field equations 

Gf,„ = df,(j) dycj) - -g^^ud^ct) d^'(j) - V{(j))g^,, + f^^ (6) 

□^-^ = -^f (7) 

where 

^-;^^ '«) 

and T = g'^'^Tn,,. Covariant differentiation of (6) and the Bianchi identities give 

V%, = f3fd,<p (9) 

which implies that the matter field is not generally conserved and feels a new force due to 
gradient of the scalar field. Let us consider T^^ as the stress-tensor of dust with energy density 
p in the Einstein frame. In a static and spherically symmetric spacetime the equation (7) gives 

df"^ f df d(j) 

where f is distance from center of the symmetry in the Einstein frame and 

V;//(0) = V{<P) - \pe-'^^ (11) 

Here we have used the relation p = e^'^^'^p that relates the energy densities in the Jordan 
and the Einstein frames. We consider a spherically symmetric body with a radius fc and a 



constant energy density pin if < fc). We also assume that the energy density outside the body 
(f > fc) is given by pout- We will denote by ipm and ipout the field values at two minima of 
the effective potential Ve//(0) inside and outside the object, respectively. They must clearly 
satisfy V^fj{ipin) = and V^j^j{ipout) = where prime indicates differentiation of Ve//(0) with 
respect to the argument. As usual, masses of small fluctuations about these minima are given 
by rriin = [VeffiVin)P and rriout = l^effiVout)]'^ which depend on ambient matter density. A 
region with large mass density corresponds to a heavy mass fleld while regions with low mass 
density corresponds to a fleld with lighter mass. In this way it is possible for the mass fleld 
to take sufficiently large values near massive objects in the Solar System scale and to hide the 
local tests. For a spherically symmetric body there is a thin-shell condition 

where $c is the Newtonian potential at f = fc- In this case, equation (10) with some appro- 
priate boundary conditions gives the field profile outside the object [19] 

= h i^out (13) 



47r Tc r 

where M^ is mass of the object. 

The function f{R) in the Jordan frame is closely related to the potential function of the 
scalar degree of freedom of the theory in the Einstein frame. Any functional form for the 
potential function corresponds to a particular class of f{R) theories. To find a viable function 
f{R) passing Solar System tests one can equivalently work with its corresponding potential 
function in the Einstein frame and put constraints on the relevant parameters via chameleon 
mechanism. Taking this as our criterion, we will use a pure exponential potential function. 
There are two reasons for this choice. Firstly, this class of potentials arises in a number of 
physical situations. In particular, there are reports that quintessence field with exponential 
potentials can produce late time acceleration [22]. Secondly, as we will show in the following 
an exponential potential function for the scalar field corresponds to a power law f{R) theory. 
In fact, there are arguments concerning cosmological viability of this class of f{R) models 
[23] [24]. Moreover, it is argued that some power law f{R) theories may have sting/M-theory 
origin [25]. These arguments make investigation of viability of these models in terms of local 
experiments be a mandate. 
To have a pure exponential potential function, we take 

where a and n are constant parameters. It then leads to 

y(0) = — i--[a !i±l]-« e2("-i)/^<^ (15) 

On the other hand, this choice of the r{p) function gives a power law f{R) theory 

f{R) = a R^+^ (16) 



Thus there is a correspondence between a power law f{R) theory in the Jordan frame and a 
minimally coupled scalar field theory with an exponential potential in the Einstein frame-'-. We 
can now find the solution of Vgjj(</)) = by substituting (15) into (11) 

1 , r n + 1 . n+ 1,„, , ^, 

09 = — nip (a )"} (17) 

In order that if he a. local minimum we should have 

n + 1 , n + L „ , , 

(a r > 18 

1 — n n 

In the following we shall consider thin-shell condition and the constraints set by equivalence 
principle and fifth force experiments. 

1. Thin — shell condition 

In the chameleon mechanism, the chameleon field is trapped inside large and massive bodies 
and its infiuence on the other bodies is only due to a thin-shell near the surface of the body. 
The criterion for this thin-shell condition is given by (12). If we combine (12) and (17) we 

obtain 

Arc 1 1 , Pout 



In ^-^ f 19) 

f, I2(3\n + l)^c Pin ^ ' 

where pm and pout are energy densities inside and outside of the body in the Jordan frame, re- 
spectively, and <l>c = Mc/87rfc with Mc being the mass of the body. In weak field approximation 
the spherically symmetric metric in the Jordan frame is given by 

ds^ = -[1 - 2X{r)\dt^ + [1 + 2Y{r)]dr'^ + r'^dU'^ (20) 

where X{r) and Y{r) are some functions of r. There is a relation between r and f so that 
f = p^/^r. If we consider this relation under the assumption niout r ^ 1, namely that the 
Compton wavelength m~^i is much larger than Solar System scales, then we have f ^ r. In 
this case, the chameleon mechanism gives for the post-Newtonian parameter 7 [27] 

3 _ Ar, 2 Ar 

7=:^^^-!-^^ (21) 



We can now apply (19) on the Earth and obtain the condition that the Earth has a thin-shell. 
To do this, we assume that the Earth is a solid sphere of radius Re = 6.4 x 10* cm and mean 
density pe ~ 10 gr /cm?. We also assume that the Earth is surrounded by an atmosphere with 
homogenous density pa ~ 10^^ grjcm? and thickness lOO/cm. Then we rewrite equation (19) 

Ai?e 1 1 , Pa 

= In — (22) 

i?e 12/32(n+l)$e Pe ^ ^ 



■'"It should be remarked that in the end of the present "work wc became a-ware of [26] in -which a relation 
bet-ween po-wer law f(R) theories and a minimally coupled scalar field -with specific exponential potential has 
been reported in the absence of matter systems. 



With $e = 6.95 X 10 ^° [28], Newtonian potential on surface of the Earth, it gives 

Ai?e _ 4.96 X 10^ 
Re ~ n + 1 



(23) 



The tightest Solar System constraint on 7 conies from Cassini tracking which gives | 7 — 1 | < 
2.3 X 10~^ [18]. This together with (21) and (23) yields an upper bound for the parameter n 

|n+l| > 1.44X 10^^ (24) 

Combining this result with (16) reveals that this power law f{R) theory hardly deviates from 
general relativity. 

2. Equivalence principle 

We now consider constraints coming from possible violation of weak equivalence principle. 
We assume that the Earth, together with its surrounding atmosphere, is an isolated body 
and neglect the effect of the other compact objects such as the Sun, the Moon and the other 
planets. Far away the Earth, matter density is modeled by a homogeneous gas with energy 
density pc ~ 10~'^'^gr/cni^. To proceed further, we first consider the condition that the atmo- 
sphere of the Earth satisfies the thin-shell condition [19]. If the atmosphere has a thin-shell 
the thickness of the shell (ARa) must be clearly smaller than that of the atmosphere itself, 
namely ARa < Ra, where Ra is the outer radius of the atmosphere. If we take thickness of 
the shell equal to that of the atmosphere itself ARa ~ 10^ km we obtain -^ < 1.5 x 10~^. 
It is then possible to relate ^^ = '^qZ^^ and ^^ = '^^Z^"' where (pe, ^a and (pc are the field 
values at the local minimum of the effective potential in the regions r < Re , Ra > r > Re and 
r > Ra respectively. Using the fact that newtonian potential inside a spherically symmetric 
object with mass density p is $ oc pR^, one can write <l>e = 10'^ $a where $e and ^a are 
Newtonian potentials on the surface of the Earth and the atmosphere, respectively. This gives 
ARe/Re ~ 10^'^ ARa/ Ra- With these results, the condition for the atmosphere to have a 
thin-shell is 

^ < 1.5 X 10-^ (25) 

Re 

The tests of equivalence principle measure the difference of free-fall acceleration of the Moon 
and the Earth towards the Sun. The constraint on the difference of the two acceleration is 
given by [18] 

\0'm ~ fle| in-13 



< 10"''' (26) 

ttN 

where a^ and Oe are acceleration of the Moon and the Earth respectively and ojv is the 
Newtonian acceleration. The Sun and the Moon are all subject to the thin-shell condition [19] 
and the field profile outside the spheres are given by (13) with replacement of corresponding 
quantities. The accelerations a^ and ag are then given by [19] 

ae^a^{l + m\^f^} (27) 



Ml + 18/3^(^f^} (28) 



-10 -T. _ Q 1 /I w 1 n-11 



where $e = 6.95 x 10"^", $„ = 3.14 x 10"^^ and ^s = 2.12 x 10"*' are Newtonian potentials 
on the surfaces of the Earth, the Moon and the Sun, respectively [28]. This gives a difference 
of free-fall acceleration 

'" ~"' (0.13) /3^ (^)2 (29) 



Combining this with (26) results in 



Rp 



< 6.74 X 10-^ (30) 



which is of the same order of the condition (25) that the atmosphere has a thin-shell. Taking 
this as the constraint coming from violation of equivalence principle, we obtain 

|n+l |> 1.67X 10^^ (31) 

which is not much different from the bound given by (24). 

3. Fifth force 

The potential energy associated with a fifth force interaction is parameterized by a Yukawa- 
type potential 

U{r) = -a— (32 

Stt r 

where mi and 1712 are masses of the two test bodies separating by distance r, a is strength 

of the interaction and A is the range. Thus fifth force experiment constrains regions of {a, A) 

parameter space. These experiments are usually carried out in a vacuum chamber in which 

the range of the interaction inside it is of the order of the size of the chamber [19], namely 

A ~ Rvac- The tightest bound on the strength of the interaction is a < 10~^ [29]. Inside the 

chamber we consider two identical bodies with uniform densities pc, radii r^ and masses rric- If 

the two bodies satisfy the thin-shell condition, their field profile outside the bodies are given 

by 

,, , l3 SArcirice ''/■^'""= ,^^, 

(Pi'^') = --A ^^vac (33) 

Air Tc r 

Then the corresponding potential energy of the interaction is 

V(r) - -2PH^? g^ (34) 

The bound on the strength of the interaction translates into 

2(5\^^f < 10-3 (35) 

One can write for each of the test bodies 

A^c _ 1 1 , Pvac f s 

r, -12/32(^ + 1)$/^ p, ^"^^^ 



where p^ac is energy density of the vacuum inside the chamber. In the experiment carried 
out in [29], one used a typical test body with mass nic ~ 40^^ and radius r^ ~ 1cm. These 
correspond to p^ ^ d.bgr/cm^ and <l>c ~ 10'^^. Moreover, the pressure in the vacuum chamber 
was reported to be 3 x 10~^ Torr which is equivalent to p^ac ~ 4.8 x 10~^^ gr/cm^^. Substituting 
these into (36) and combining the result with (35) gives the bound 

|n + l|>3xl02^ (37) 

which is much stronger than (24) and (31). As the last point, there are some remarks to do 
with respect to stability of the model (16). In principle, stability issues should be considered 
to make sure that an f{R) model is viable [30]. In particular, stability in matter sector (the 
Dolgov-Kawasaki instability [31]) imposes some conditions on the functional form of f{R) 
models. The first theories which easily pass this instability have been presented in [20] [32]. 
These conditions require that the first and the second derivatives of f{R) function with respect 
to the Ricci scalar R should be positive definite. The positivity of the first derivative ensures 
that the scalar degree of freedom is not tachyonic and positivity of the second derivative tells 
us that graviton is not a ghost. For power law f{R) theories of the type (16), we should have 
n > —1 and n > to ensure that f'{R) > and f"{R) > 0, respectively. 
In summary, we have discussed viability of power law f{R) theories in terms of local gravity 
constraints. We have used the correspondence between this class of f{R) theories to scalar 
field theories with an exponential self-interacting potential. In the scalar field representation 
of the theory there is a strong coupling of the scalar field with the matter sector. We have 
considered the conditions that this coupling is suppressed by chameleon mechanism. We have 
found that in order that the theory be consistent with local gravity experiments the exponent 
of the curvature scalar hardly deviates from unity. The constraint (37) is much stronger than 
that reported in [33] which is obtained by considering the perihelion precession of Mercury 
under the assumption that it follows timelike geodesies. Our results preclude the possibility of 
regarding power law f{R) models as viable candidates for generalizing general relativity. 



^Note that due to the logarithmic dependence of -^ with p^ac the value of the vacuum energy density does 
not effectively change the order of magnitude of (37). 
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